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Forced-oscillation  pitch-damping  balances  oscillate  over  a 
small  angular  amplitude  range  about  a  nominal  angle  of  attack  -and-thus 
yield  an  effective  value  of  the  aerodynamic  damping  if  the  damping  is 
a  nonlinear  function  of  angle  of  attack.  Because  the  local  value  of 
the  damping  coefficient  is  generally  required  for  computer-simulated 
motion  studies,  a  procedure  to  extract  the  local  damping  from  the 
experimental  effective  damping  output  of  a  forced-oscillation  balance 
is  derived.  A  derivation  is  given  of  the  basic  integral  equation 
relating  local  and  effective  damping.  Techniques  for  solving  this 
integral  equation  are  given.  The  method  is  applied  to  experimental 
nonlinear  damping  data  for  three  preliminary  unmanned  Mars  reentry 
(Viking)  configurations,  and  the  results  are  discussed. 
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ABSTRACT 


Forced-oscillation  pitch-damping  balances  oscillate  over  a  small  angular  amplitude 
range  about  a  nominal  angle  of  attack  and  thus  yield  an  effective  value  of  the  aerodynamic 
damping  if  the  damping  is  a  nonlinear  function  of  angle  of  attack  Because  the  local  value 
of  the  damping  coefficient  is  generally  required  for  computer-simulated  motion  studies 
a  procedure  to  extract  the  local  damping  from  the  experimental  effective  damping  output 
of  a  forced-oscillation  balance  is  derived.  A  derivation  is  given  of  the  basic  integral  equation 
rcla'mg  local  and  effective  damping.  Techniques  for  solving  this  integral  equation  are  given 
The  method  is  applied  to  experimental  nonlinear  damping  data  for  three  preliminary 
unmanned  Mars  reentry  tViking)  configurations,  and  the  results  are  discussed. 
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NOMENCLATURE 

Ak  k,h  Fourier  coefficient  for  Cm  D,  Eq.  (33)  and  Fig.  5 

ak  k,h  Fourier  coefficient  for  CmDt,f,  Eq.  (33)  and  Fig  5 

B„  Balance  flexure  spring  constant,  ft-lb/rad 

B.'t  Balance  damping  tare,  ft-lb-scc/rad 

b  Slope  defined  by  Eq.  (61),  per  rad 

C0-C3  Coefficients  in  the  polynomial  representation  of  Cm  D  and  Cm  D  (Eqs  (34) 
and  (39)) 

Cm  Pitching-moment  coefficient,  pitching  momcnt/q^Sd 

On  £)  Cm  q  +  Cm  ’a ,  per  rad 
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fma  3Cm /3a,  per  rud 
C,„;  3Cm/3(ad/2V.),  par  rud 

Cmq  3Cm /3(qd/2V„),  per  rad 

d  Maximum  medal  diameter,  ft 

EVENS  Defined  by  Eq.  (35) 

I  Model  moment  of  inertia  about  pivot  axis,  ft-lb-scc/rad 

t  j,k  General  subscripts  and  indices 

12  Modal  length,  ft 

Mj  Total  aerodynamic  moment,  Eq  (20),  ft-lb 

Mi  Defined  by  Eq.  (17),  ft-lb 

Mo  Mq  +  Mi ,  ft-lb-sec/rac! 

M(t)  Forcing  moment  of  balance,  ft-lb 

Mj,  Hysteresis  moment,  ft-lb  ( Eq  (28)) 

M_  Free-stream  Mach  number 

M„  Aerodynamic  pitching-moment  slope,  q_Sd  C„,a.  ft-lb/rad 

Mi  ((q„Sd2)/(2V„)|  Cn;  .  ft-lb-scc/rad 

M,  [(q_Sd2)/(2V_)j  Cmq.  ft-lb-scc/rad 

My  Sum  of  aerodynamic,  balance  flexure,  and  forcing  moments  acting  on  model, 

ft-lb  (Eqs  (l)’and  (2)) 

N  Indices  for  a  where  Cm[)  =  Cmi)  (Fig  4) 

ODDS  Defined  by  Eq  (36) 

Pj-P8  Defined  by  lap  (49)  through  (56).  respectively 

q  Pitching  rate,  in  general,  and  a.  in  particular,  for  the  balance,  rad/sac 

o.  Free  st.,  m  dynamic  pressure.  Ib/it2 
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Ro-Rj  Defined  by  Eqs.  (40)  through  (48),  respectively 
Ree  Frce-stream  Reynolds  number  based  on  model  maximum  diameter 

S  Reference  area,  jrd2/4,  ft2 

t  Time,  sec 

Tt  k'h  coefficient  in  series  representation  of  M(t)  ft-Ib 

V_  Free-stream  velocity,  ft/sec 

Wk  k<h  weighting  function  relating  Ak  to  ak ,  Eq.  (33) 

x<t  Model  center  of  gravity  location  with  respect  to  nose,  in.  (Fig.  9) 

a  Instantaneous  angle  of  attack  of  model,  rad 

ag  Defined  in  Fig.  8,  rad 

amix  a  limit  in  Fourier  series  representation  of  CmD ,  rad  (Fig.  S) 
aN  a  where  Cm  „  =  Cm  D<;(f  (see  Fig.  4),  rad 

do  Nominal  or  equilibrium  angle  of  attack,  Eq.  (4),  rad 

aST  Balance  sting  angle  of  attack,  rad 

a  Time  rate  of  change  of  angle  of  attack,  rad/sec 

0  Integral  limit  defined  by  Eq  (37),  rad 

7t.72.73.  Integral  limits  defined  by  Eqs.  (57),  (58),  and  (59),  respectively,  rad  (Fig.  8) 
Aa  Increment  of  a  defined  in  Fig  4,  rad 

0  Instantaneous  angular  oscillation  position,  rad 

0  Time  rate  of  change  of  6,  rad/sec 

Oj  Maximum  angular  oscillation  amplitude  for  a  cycle  of  motion,  rad 


8o,6„0c  Defined  in  Fig.  2,  rad 

P  Phase  angle  between  forcing  moment.  M(t),  and  oscillation  amplitude.  8(t),  rad 
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of  M(t)  or  0(t),  rad/sec 
of  balance  system  V7m„  +  Ua  )l\,  rad/sec 

SUBSCRIPTS 

eff  Effective  value 

ST  Sting 


Angular  circular  frequency 
Natural  circular  frequency 
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SECTION  I 
INTRODUCTION 

Two  general  procedures  are  available  for  measuring  pitch  damping  dynamic  stability 
coefficients  for  a  captive  (sting-supported)  model  in  a  wind  tunnel  (Ref  !)  These  methods 
are  the  forced-oscillation  and  the  frce-osciilation  techniques  The  frec-oscillation  balance 
is  usually  somewhat  simpler  n  design  and  operation  than  a  forced-oscillation  balance,  but 
the  forced-oscillation  balance  provides  more  precise  control  of  the  amplitude  of  the 
oscillation,  which  can  be  particularly  important  if  the  model  is  dynamically  unstable  Data 
reduction  for  a  forced-oscillation  system  is  generally  simpler  than  for  a  frec-oscillation 
balance.  The  free-ose illation  balance  yields  transient  angular  motion  data  which  must  be 
analyzed  in  terms  of  the  decay  in  the  amplitude  to  yield  the  desired  damping  coefficients 
As  described  below,  the  forced-oscillation  balance  operate*,  at  constant  angular  motion, 
and  the  damping  coefficients  are  obtained  directly  troni  the  torque  input  measurements 

The  von  Karmen  Gas  Dynamics  Facility  (VKF)  of  the  Arnold  Lngmeenng 
Development  Center  (ALDC)  has  developed  a  forced-oscillation  dynamic  balance  to  measure 
pitch  (or  yaw)  damping.  This  balance  is  sting  supported  with  a  cross-flexure  spring  pivot 
system  which  allows  one  degree  of  freedom  in  angular  motion  Pie  balance  system  is 
described  in  Refs.  1  and  2.  Certain  essential  features  of  the  system  operation  are  also 
described  herein  to  provide  the  necessary  background  for  the  present  analysis  A  schematic 
representation  of  the  VKF  forced-oscillation  balance  is  depicted  in  Fig  1 

The  balance  is  equipped  with  a  magnetic  shaker  motor  and  feedback  control  (Fig 
2)  to  provide  oscillation  amplitudes,  0j,  from  essentially  zero  to  approximately  2  deg 
However,  to  reduce  effects  of  tunnel  vibrations,  it  is  generally  necessary  to  make  the 
oscillation  amplitude  greater  than  about  1  deg  The  balance  system  output  is  a  damping 
coefficient  which  is  the  mean  value  over  a  cycle  of  motion  for  that  particular  oscillation 
amplitude,  0 j .  at  the  nominal  angle  of  attack.  a0.  As  the  oscillation  amplitude  is  lowered, 
the  "effective"  damping  determined  by  the  system  approaches  the  "local"  value  for  that 
particular  q^.  If  the  damping  coefficient  is  a  linear  function  of  a  within  the  region  a,, 
+  01  to  Oq  -  0i ,  then  the  balance  system  gives  the  true  local  value  of  the  damping  coefficient 
on  that  particular  do. 

The  values  of  the  damping  coefficient  obtained  from  the  forced-oscillation  technique 
arc  easily  interpreted  if  the  variation  of  the  local  value  of  the  damping  coefficient  with 
angle  of  attack  is  small  within  the  range  covered  during  the  oscillation  from  a0  -  0i 
to  do  +0i.  There  exist  cases  for  which  the  variation  of  damping  coefficient  over  the 
range  of  angle  of  attack  is  quite  large.  One  such  case  was  the  unmanned  Mars  landing 
vehicle  (Viking),  tested  in  the  AHDC  Propulsion  Wind  Tunnels  (Refs  3  and  4),  which 
was  found  to  be  dynamically  unstable  for  small  angles  of  attack  (less  than  about  2  deg) 
and  dynamically  stable  at  higher  angles.  The  variation  of  the  damping  coefficient  with 
angle  of  attack  over  the  angle-of-attack  range  covered  in  the  forced  oscilljtion  was  quite 
large.  The  purpose  of  tins  report  is  to  describe  the  mathematical  techniques  by  which 
the  value  of  the  damping  coefficients  for  very  small  oscillations  (termed  the  local  value, 
Cm  ()  (do))  can  be  determined  from  the  effective  value  obtained  with  the  forced  oscillation 
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Fig.  1  Seth.  iatic  Representation  of  Balance  System 
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„  „  (Static  Moment 
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Bj  (tt  ■  Error  Signal  •  Bj  (t)  -  Bg  It) 

Fig.  2  Schematic  of  Balance  Feedback  Control  System 
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balance  by  oscillating  the  system  with  an  amplitude,  0j.  about  the  angle  of  attack  Oo, 
Cm  D  (ao.O,)  As  the  value  of  Q\  goes  to  zero,  the  effective  value  and  the  local  value 
arc  the  same;  i  e., 

0,lim  Cm  («.0,)  -  Cm  (a) 

0  mU  „  °  1  ml)  ° 

fii 

The  local  damping  coefficients  arc  of  practical  importance  because  they  arc  required  in 
computer  simulations  of  dynamic  motions. 

SECTION  II 
ANALYSIS 

2.1  FORCED-OSCILLATION  BALANCE  SYSTEM  OUTPUT 

A  schematic  of  the  VKF  forced-oscillation  balance  is  shown  in  Fig.  1  A  magnetic 
shake,  motor  is  used  to  impose  an  approximately  sinusoidal  oscillation  upon  the  model 
The  moments  opposing  the  moment  input  by  the  motor  are  those  attributable  to  inertia, 
the  aerodynamic  restoring  moment,  and  the  flexure.  Strain  gages  arc  attached  to  the  input 
toiquc  beam  and  the  flexure.  The  first  of  these  gives  the  input  momenj  from  the  shaker 
motor,  and  the  second  is  used  to  determine  the  angular  position  of  the  model.  0. 


A  feedback  into  the  shaker  motor  control  is  used  to  maintain  the  amplitude  of  the 
model  oscillation.  0|,  at  a  preselected  value  (Fig.  2).  It  is  this  feedback  that  allows  the 
device  to  be  used  even  with  dynamically  unstable  systems  since,  for  these,  the  shaker  motor 
extracts  energy  from  the  system  being  tested  With  a  dynamically  stable  system,  the  shaker 
motor  must  add  energy  to  the  system  to  maintain  a  constant  ampi.'udc. 

To  determine  the  damping,  the  frequency  is  varied  until  a  phase  shift.  <y  of 
approximately  90  deg  between  the  input  torque  and  the  displacement.  P,  is  obtained 
At  this  condition  the  moments  attributable  to  inertia  and  stiffness  arc  essentially  in  balance, 
so  the  output  of  the  shaker  motor  is  basically  compensating  for  the  damping  moment. 
Tlic  amplitudes  of  the  moment.  Tj.  and  the  displacement,  0i .  along  with  the  frequency, 
w.  allow  the  damping  term  to  be  evaluated 


.m 


The  control  console  is  shown  schematically  in  the  lower  portion  of  Fig.  1.  Internally, 
there  arc  low  pass  filters,  phase  shifters,  and  multipliers,  winch  arc  used  to  give  the 
amplitudes  of  the  in-phase  and  out-of-phase  components  of  the  input  torque  From  these 
the  aerodynamic  damping  and  stiffness  may  be  determined. 

The  basic  differential  equation  for  a  forced-oscillation  pitch-damping  balance  is 


Mjja  *  Mfla  *  B,j«i  * 


Arrochnjmw 
Hectoring  Moment 


H.dorur  1  !r\urr 
MomrPl 


1  orctng  Moment  (-) 

of  ''holer 
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I  =  Model  moment  of  inertia  about  pivot  axis 

Mr  =  Sum  of  aerodynamic,  balance  flexure,  and  forcing 
moments  acting  on  model  about  pivot  axis 

Mp  =  Aerodynamic  damping  moment  per  angular  velocity 
=  (Cm;  +  C„q)  q.SdV2V. 

M0  =  Aerodynamic  pitching  moment  per  unit  angle  (stiffness) 
=  Cn„  fl-Sd 


=  Balance  flexure  structural  damping  moment  per 
angular  velocity 

B„  =  Balance  flexure  moment  per  unit  angle  (stiffness) 


M(t)  =  Forced  oscillatory  moment  impressed  on  model 
asT  =  Stationary  angle  of  attack  of  sting 

a  =  Instantaneous  angle  of  attack  of  the  model 


The  angle  of  attack  may  be  written  as 


a  «  a  ♦  d(l) 


0(t)  =  Instantaneous  oscillation  amplitude  about  Og 
a,,  =  Equilibrium  angle  of  attack  of  mode! 


The  relation  between  the  sting  angle  and  the  equilibrium  angle  is  given  by 


vVST 


The  relationships  among  the  various  angles  are  shown  in  Fig.  3. 

After  substituting  Eqs.  (2),  (3),  and  (4)  into  Eq  (1),  the  basic  balance  equation 
becomes 


AEOC-VR72-25 


-  Sting  Centerline 


-  Model-Balance  Centerline  for  Static 
Air-on-Equilibrium  Conditions 


Maximum 
Oscillation 
^  Amplitudes 


~T 

Relative  Wind  Direction 
Balance  Flexure  Apparent  Center  of 
Rotation  and  Model  Center  of  Gravity- 


Fig.  3  Geometry  of  Pitch-Damping  Balance-Sting  Combination 


The  stiffness  of  the  flexure,  B„,  is  generally  selected  to  be  much  greater  than  the 
aerodynamic  stiffness,  M„,  since  the  latter  varies  widely  with  tunnel  conditions.  This  also 
decreases  the  effects  of  nonlineantics  in  the  static  restoring  moment  In  addition,  the 
system  is  designed  so  that  Ba  is  constant  over  the  range  of  deflections  encountered  This, 
in  conjunction  with  the  feedback  control  system,  permits  the  oscillation  amplitude  to 
be  accurately  represented  by  a  pure  cosine  wave 

fl(i)  =  t?|  (5) 

Tne  output  of  the  torque  beam  will,  in  general,  contain  higher  harmonics,  since  by 
Eq.  (5)  ail  the  damping  and  nonlinear  terms  will  be  included  in  M(t).  Thus,  the  input 
moment  can  be  expressed  in  a  Fourier  scries 

Mft)  «  io  -  <\i,|.  lyo  -  'm(  ,  ::a)o 

f  7) 

w  ^  co*  (Ul  *  <5) 

Ul 


If  all  of  the  coefficients  in  the  basic  differential  equation  are  constants,  then  only 
the  fust  term  (k  =  I)  is  present  in  M(t).  Using  Eq  (6)  and  collecting  the  coefficients 
of  cos  tot  (in  phase)  and  sin  cot  (out  of  phase)  give  for  this  particular  case 


-  (Ma  ,  \\(l)0i  -  T,  co 

-  (Mjj  ♦  ^  \  s,n  ^ 


(8) 


5 


JK** 
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Thus 


51  ♦  B, 


— hu^  —  ^-2  cos  0 


equals  total  static  moment  slope,  and 


'<D 


equals  total  system  damping  coefficient. 


(9) 


(10) 


Normally  the  frequency  is  selected  (<o  =  <o„)  to  give  a  phase  shift  of  90  deg  in  order 
to  maximize  the  accuracy  with  which  the  damping  coefficient  may  be  determined. 

...  Tlle  rncasure<1  values  of  T( ,  *,  to,  and  0\ ,  and  values  of  B,  and  B;  from  pretest 
ralihrations  then  aHow  the  aerodynamic  coefficients,  M„  and  MD,  to  be  determined  from 
tqs.  (9)  and  (10).  Normally,  the  aerodynamic  stiffness,  M„,  may  be  determined  more 
accurately  with  a  static  balance.  The  aerodynamic  damping,  MD ,  is  the  sum  of  M  •  and  M„ . 
The  individual  contributions  of  and  M,  cannot  be  distinguished  by  the  balance,  since  a 


2.2  EQUIVALENT  DAMPING  FOR  NONLINEAR  SYSTEMS 

If  there  are  nonlinearities  in  the  system,  either  or  both  of  the  coefficients  MD  and 
may  vary.  In  general,  each  may  be  a  function  of  a  and  a.  To  avoid  certain  pitfalls, 
it  is  preferable  to  begin  with  the  expression  for  the  total  aerodynamic  moment  M,  as 
a  function  of  a  and  a 


"a  *  X,a(a*«)  (II) 

The  circuitry  of  the  console  contains  low-pass  filters  so  that  the  higher  harmonics 
of  M(t)  aie  filtered  out.  Therefore,  the  system  gives  the  fundamental  harmonic;  i.e.,  the 
m-phasc  (Tj  cos  0)  and  the  out-of-phase  (T1  sm  0)  fundamental  components  of  the  input 
torque.  In  genera! 


and 


M(t)  -  10  -  Bd0-  BQ0  -  M Jao  +  et0) 

■  £  cos  (kwt  +  0) 
k.l 


T|  COS  “  J M(|)  coscot  d(o)t) 
o 

2  n 

-T|  sin  0  .  ^M(t)  sin  cut  cKcvl) 


(12) 

(13) 

(14) 


6 


mammam 
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With  0  ~~  d\  cos  cot  these  become 


1|COS^  =  '"fl-lu-O,  cosut  ♦  B^fl,  s.nn>t  _  iia0)CO5W|  -  M  j  COS  cut  (Knit)  05) 
1  r2n 

—  1 1  sin  d  ,  -J(-|w=9|cosO,,Bio<?l8.n<u.-Ba0!co8«Jt_Jlj,I„4),<i(o,)  (16) 

with 


^  =  '!>„  ♦  0|  eo*  oil,  —  oiO j  sin  cut) 

Certain  of  the  integrals  may  be  evaluated,  with  the  result 


(17) 


Tl  COS  6  -  -hrOl  -  | i  0t 


08) 


-T|  sin  <J  »  B^-oifl, 


<  9,  coswt ,~oOl  sin  ot)  sin  mi  cKcjt) 


09) 


These  results  are  valid  for  arbitrary  aerodynamic  moments. 


mnmnnf ^'0f  '°  STpli,V  thcse  rcsults  for  certain  specific  cases  If  tile  aerodynamic 

moment  is  a  continuous  function  of  the  angular  velocity,  then  an  expansion  can  be  made 
about  an  angular  velocity  of  zero. 


M„(a,a) 


<?\l  dim 

'<>.0>  *  ~ia'  0)i  *  ’* - -  (a  ,  Old 2  *  ... 

da  d*i 


In  terms  of  these  quantities 


\  -  M>,0Va 


M„ 


(V  (a.a)-M  (a. 0)1 
• _ * _ * 

a 


9M.  9SM 

*  fa.O)  -♦  !j  — fa,0)a  -»  ... 
di  di¬ 


tto) 


(2i) 


(22) 
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The  parameter  M„  is  therefore  independent  of  velocity.  This  is  a  direct  consequence  of 
the  assumption  that  Ma  is  continuous,  since  any  dependence  upon  velocity  appears  in 
Mu .  The  moment  can  r.ow  be  written  as 


M  Jata)  -  Ma(ala  +  MQ{a,d)d 


With  this  representation,  the  integrals  may  be  simplified 

T,  cos  ^ 

- 2 - m  -Iw  —  B  —  — r—  J  M  (a)a  coscjt  cKc^t) 

V  J  a  m) ,  ▼  a 

f'n 

■  "Iw  “  “”]  Ma(<*o  ^  cos  ^  C08^  <Ka>0 


(23) 


(24) 


I  .  MQ  O  \e  *> 

-~g-  ■  »«♦  Tj  cos o)t,  —  cj6 j  sin  <j0  sin^ctil  <Xa>t)  (25) 

*  o 

According  to  this,  nonlinearities  in  the  static  moment  do  not  affect  the  damping  as 
determined  by  the  forced-oscillation  balance.  This  is  valid  if  Ms  is  a  single  valued  function 
of  a. 


If  the  damping  coefficient  is  expanded  in  terms  of  a,  i.e. 


M  |)(a,a) 


MD(a,0) 


- (a,0)a  ♦ 

da 


(26) 


Equation  (25)  becomes 


T ,  *in  <}>  j  •  2rr  ^ 

- ^  ~  J  ^l|)(ao  ♦  d|  cos  cjt,  0)  sin*  o>t  <Kwt) 


,217  d2M„ 

•  “  j  - tag  +  (?j  cos6>t,  0)  sin4  cjt  (27) 

•  da 

to40j  *,2n  d4MD 

+  — —  i  —— 7-  (aft  +  0.  cosot,  0)  sin**  o>t  d(cut)  *■ 

17  da*  °  1 

The  left-hand  side  of  this  equation  is  the  value  determined  by  the  forced-oscillation 
technique.  If  the  aerodynamic  damping  coefficient  Mp  is  constant  or  is  a  linear  function 
of  a,  then  the  result  reverts  to  that  given  by  Eq.  (10).  If  the  aerodynamic  coefficient 
is  nonlinear  with  angular  velocity,  i.e. 
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according  to  Eq.  (27)  the  apparent  effective  damping  determined  by  the  forced-oscillation 
technique  will  depend  upon  the  frequency,  w.  For  this  reason  it  is  advisable  to  obtain 
data  at  two  or  more  frequencies,  which  requires  using  more  than  one  set  of  flexures 
If  there  is  no  dependence  upon  frequency  noted,  then  it  may  be  concluded  that 


d2nv 


D 


da 


2n 


0 


One  class  of  nonlineanties  that  does  not  fall  within  the  above  framework  is  hysteresis, 
because  the  moment  is  not  a  continuous  analytic  function  of  velocity.  If  Mh  is  the  hysteresis 
moment,  and  assuming  for  simplicity  that  the  variation  of  Ma  with  a  is  linear 


M,  -  M„a  ♦  Mk  a  >  0 

•i  a  n 

"  '!„«  -  "h  «  <  0 


(23) 


Then 


T|  sin  0 

(dO, 


’ll 

o>0. 


(29) 


In  this  case  the  apparent  damping  varies  inversely  with  udt.  Since  the  balances  operate 
at  constant  c o0\ ,  the  quantity  (4Mh/jra)0|  )would  appear  to  the  balance  system  as  a  constant 
value  of  effective  aerodynamic  damping,  Mpc(f.  for  that  particular  test  condition.  As  with 
the  case  of  a  nonlinear  damping  coefficient,  testing  at  different  frequencies  would  rcvcai 
the  existence  of  the  effect. 


In  the  remaining  development  it  is  assumed  that  there  is  no  hysteresis,  and  the 
nonlineanties  with  respect  to  angular  velocity  in  the  damping  coefficient  arc  negligible 
Under  these  conditions 

t  A  i  r1" 

- —  =  11^  ♦  —J  Mj)(a0  *  0 1  costal,  0)  sin ‘tat  d(tal)  (30) 

*  o 

If  now  the  local  value  of  damping  derivative,  Cm  D  (a),  is  related  to  Mp  and  the  effective 
value,  C„D[|[  (Oo,  8| ),  is  that  determined  from  the  forced-oscillation  technique,  i  c  ,  from 
tq.  (30),  then 


If2”  . 

C,  (0o.d,)  -  —  J  Cy  (an  +  0,  cos  cal)  sin*  tat  dt 
.ft  t>  ** 


(31) 


Equation  (31)  is  an  integral  equation  which  relates  the  effective  cxpcnmcntal  damping 
obtained  from  the  present  balance  system  to  the  local  damping  Tins  integral  relationship 
can  also  be  obtained  directly,  without  considering  the  balance  output,  by  equating  the 
amount  of  energy  dissipated  per  cycle  by  the  local  and  effective  damping  (Refs  5,  6. 
and  7). 
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Equation  (31)  is  valid  for  a  nonlinear  pitching  moment  (Cm  versus  a)  so  long  as 
the  pitching  moment  is  single  valued  for  a  given  a.  That  is,  pitching-moment  hysteresis 
or  time  lags  (because  of  viscous  effects)  can  also  add  or  dissipate  system  energy,  which 
will  then  appear  as  nonlinear  damping  (Eq  (29)). 

23  SOLUTION  OF  INTEGRAL  EQUATION 

Equation  (31)  is  a  Fredholm  integral  equation  of  the  first  kind  (Ref.  8).  The  known 
quantity  is  CmBeff  which  is  a  function  of  and  6\  and  usually  is  known  for  various 
values  of  Oo  at  a  fixed  value  of  0|.  The  unknown  quantity  is  CmD,  which  is  a  function 
of  a.  Three  ways  of  solving  this  equation  were  tned,  with  the  first  two  being  unsuccessful. 

For  the  first  approach,  Eqs.  (3)  and  (6)  were  used  in  Eq  (31)  to  give 

■  A  <3S 


Equation  (32)  was  written  in  finite  difference  form  involving  the  known  values  of  Cn()c[f 
and  the  unknown  values  of  Cm  0  (Fig.  4)  The  procedure  was  very  similar  to  that  mentioned 
in  Ref.  9.  This  formulation  led  to  a  set  of  linear  simultaneous  equations  for  the  unknown 
quantities  which  could  be  solved  by  standard  techniques.  The  answers,  however,  were 
extremely  sensitive  to  the  input  values  of  the  effective  damping.  The  results,  in  general, 
were  very  erratic,  even  for  simple  test  cases,  and  the  method  was  abandoned. 

The  second  approach  to  obtain  a  practical  relationship  between  the  local  and  effective 
damping  consisted  of  a  scries  solution  for  Eq.  (32)  (Ref.  8)  The  known  effective  damping 
and  the  unknown  local  damping  were  represented  by  Fourier  cosine  senes.  The  kernel 
was  represented  by  various  types  of  power  series  in  a.  This  allowed  the  analytic  evaluation 
of  the  integral  for  a  general  term,  and  a  relation  between  the  known  (effective)  and 
unknown  (local)  Fourier  coefficients  could  be  established.  This  relation  for  the  k,h 
coefficient  is  of  the  form 


The  weighting  coefficient,  Wk,  is  shown  schematically  in  Fig.  S  as  a  function  of  k.  Note 
that  for  large  values  of  k,  Wk  is  very  small  and  oscillates  about  zero.  It  was  this  osciliatory 
behavior  of  Wk  that  led  to  the  abandonment  of  this  method,  since  Eq.  (33)  could  become 
singular  unless  ak  was  determined  from  "mathematically  perfect"  data  for  Cm  p This 
basic  procedure  has  been  proposed  (but  not  applied)  recently  in  Ref  10,  where  Wk  is 
evaluated  by  Bessel  functions.  Unfortunately,  the  method  places  such  a  premium  on  the 
quality  of  the  input  data  and  the  Fourier  senes  representation  of  this  data  that  it  is  not 
practical  for  this  particular  case. 
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j  ■  1  -  *  N ,  a0  ■  j  A  a  ,  8j  •  k  A  a 

This  provides  N  linear  equations  which  may  be  solved  (or  N  unknown  values  o(  CmQ . 


Fig.  4  Schematic  Illustration  of  the  Damping  in  Discrete  Variable  Form 


!  I 
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Fig.  5  Typical  Variation  of  Wj,  versus  k 
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>2|K':  |  Both  of  the  above  techniques  failed  because  of  a  very  strong  sensitivity  of  tilt  local 

wH|  |  damping  coefficient  to  the  effective  damping  coefficient.  Many  physically  unrealistic  values 

jflHp  I .  of  the  local  damping  coefficient  were  obtained,  such  as  values  that  were  alternately  large 

I*  positive  and  large  negative.  The  reason  for  this  behavior  can  be  seen  in  Fig  5,  which 

|  essentially  states  that  there  are  certain  specific  values  of  higher  frequency  contribution 

|.  which  should  not  be  present  in  the  effective  damping.  Because  of  noise  and  other 

|  inaccuracies,  such  frequencies  will  be  present,  and,  in  the  inversion  process,  effects  of 

f  the  noise  will  dominate. 

i  ■ 

I-'  ' 

t 

i  'f 

l . 

!■ 

l  To  overcome  this  undesirable  effect,  smoothing  of  the  input  data  is  necessary  In 

*  the  second  approach  smoothing  by  Fourier  senes  was  found  to  be  unacceptable  The  hird 

jMt.  approach  made  use  of  polynomial  smoothing.  The  manner  in  which  this  was  accompli  died 

*'  is  somewhat  similar  to  the  procedure  given  in  Ref.  6  for  free-osc illation  pitch-damping 

data,  However,  certain  features  of  the  present  analysis  for  forced  oscillation  were  not 
necessary  in  Ref.  6  and,  as  such,  ment  inclusion  here.  To  provide  continuity,  the  pro<  edure 
wBrav  is  outlined  ir.  detail 

1  .  ’. 

■f 

't  . 

It  is  assumed  that  CmD  can  be  represented  by  a  powci  scries  in  a  of  the  form 

Hf  *  (:M„  -  Co  *  C,|«|  ♦  Cia!  .  C3la|J  .  C404  .  C.  ifl|5  .  c6««  •  C.  V'  .  (.8a8  (.14) 

(MB'  The  odd  powers  are  included  to  enhance  the  curve  fitting  capability,  but  they  must  be 

considered  as  absolute  quantities  because  of  the  requirement  that  Cml)  be  symmetrica! 
ioBe  about  an  angle  of  attack  of  2ero.  This  is  valid  for  a  body  of  revolution 

^BEI  »  liquations  (3)  and  (6)  arc  substituted  into  Eq  (34).  The  resulting  expression  for 

Cm  u  may  be  substituted  into  Eq  (30)  and  integrated  analytically.  It  is  convenient  to 
,  work  with  tlie  even  and  odd  terms  separately  The  even  terms  of  Eq  (34)  arc 

|B|  E\ ENS  *  c(t  -  cosoif  *  tow4 

4  C6(«o  *  0  j  coswl)®  ♦  Cg(«o  ♦  0\ 

1  ■* 

l 

\ 

i  ■ 

HK)  The  odd  terms  of  Eq.  (34)  arc 

j  ODDS  -  C |(n0  .  dj  cos  wt)  -  C,(«„  .fl,  ...scat)1 

NM|  j  .  C-(aM  *  4|  cos  col)  ’  ‘  (..(nQ  ,  0 1  cos  coll' 

j9^B|  |’  If  eq,  is  less  than  6 1 ,  then  a  will  be  negative  over  a  portion  of  the  half  cycle  (Figs 

|  6  and  7).  Again  referring  to  Figs.  6  and  7  and  using  Eqs.  (3)  and  (6).  0  can  be  determined 

L 
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Fig.  7  Typical  Cyclic  Variation  of  a  and  C, 
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The  odd  terms  are  made  positive  (or  absolute)  in  the  cyclic  region  from  ir  -  0  to  it  by 
changing  sign.  With  rcfcience  to  Fig.  7  and  Eq.  (31),  the  same  amounts  of  energy  are 
dissipated  over  each  half  cycle  Consequently,  Eqs.  (35),  (36),  and  the  results  of  Eq.  (37) 
can  be  combined  with  Eq.  (31)  to  give 

IT 

C  *  —  f  (EVENS!  sin2  (fj t)  d(wt) 

D,.l 

*  J  [ODDS!  sin2  fnit)  (Kuli  (38) 

o 

-  |  /  [ODDS}  Sin*  (wt)  d<CO*> 

The  evaluation  of  Eq.  '38)  yields 
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»>3  " 

(SI) 

P  4  * 

Id  s'atP- 1  *i"s-s) 

(52) 

Ps  * 

2  (7T  _  hid  2/3  4/3  sin  6/3\ 

tt\2  ~  4  *  4  +  12/ 

(53) 

P6  ' 

2  /mo3  ^  2  5  Q  h,o7  /3\ 

nii-l  j-j 

(54) 

Pt  * 
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If  Ho  =  0,  the  above  results  also  apply  to  the  type  of  free  oscillation  about  zero 
a  which  was  considered  in  Ref.  6  In  particular,  if  the  odd  terms  arc  neglected  for  this 
condition,  Eq.  (39)  reduces  to  Eq.  (13)  of  Ref.  6.  The  polynomial  method  may  a1  so 
be  applied  to  frce-oscillation  test  data  for  do  not  equal  to  zero  if  the  frcc-oscillation 
damping  is  determined  for  the  same  value  of  0,  from  oscillation  envelopes  at  different 
values  of  do- 


Equation  (34)  (Cra„)  and  Eq.  (39)  (CroDcff)  are  related  in  that  both  contain  the 
coefficients  C0...C6  as  unknown  quantities  Practical  use  is  made  of  this  relationship  by 


curve-fitting  the  experimental  values  of  CmDe[f  as  functions  of  o„  and  0\  with  !.q  (.39) 


and  determining  the  coefficients  Co...C8  by  a  least-squares  procedure.  This  immediately 
yields  the  corresponding  values  of  Cn  n  since  they  can  now  be  computed  from  Eq  (34) 
for  a  given  a  =  (to. 


This  analysis  was  made  for  an  eighth-order  polynomial  curve  (nine  unknowns),  but 
it  is  easily  reduced  to  a  sixth-order  (seven  unknowns)  or  a  fourth-order  (five  unknowns) 
curve  by  merely  deleting  the  higher-order  terms.  A  program  has  been  written  for  the  VKE 
C DC-1 604 B  computer  which  incorporates  the  above  analysis  relating  the  local  and  effective 
damping  coefficients.  The  coefficients  arc  determined  for  a  minimum  square-error  fit  to 
the  experimental  data  points. 


2.4  THE  V-WELL  AND  SQUARE  WELL 


One  of  the  earlier  approaches  to  acquire  insight  concerning  the  local-effective  damping 
rclationsliip  expressed  by  Eq.  (31)  was  to  assume  simple  variat.ons  for  the  local  damping 
and  compute  the  resulting  effective  damping.  One  such  variation  was  the  V-wcll 
representation  of  the  local  damping  as  depicted  in  Fig.  8.  This  simple  variation  is  composed 
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K  s/2  0 


a.  Schematic  Depiction  of  Integral  Limits 


b.  Typical  Effect  of  Oscillation  Amplitude  on  Effective  Damping 
Fig.  8  V-Well  and  Square  Wall  Representation  of  Local  Damping 
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(62) 


c,Ky  -  =»] 

*  (t  — 4  £2 si"3  y2  ~  B'n3  y>  -  ’'"3  >3]} 

The  final  result  for  the  square  well  is: 

S,  (ao’  ■  <v^  *  l[cmD(«K>-  [(?  -  -  (7  -  -^)]  (63) 


In  a  like  manner,  other  variations  of  the  local  damping  composed  of  straight-line  (linear) 
segments  can  be  handled.  However,  the  above  results  are  sufficient  to  represent  large 
gradients  or  discontinuities  in  the  local  damping  near  zero  angle  of  attack  and  predict 
the  balance  output  under  such  conditions.  For  example.  Fig.  8b  illustrates  the  rather  large 
effect  of  oscillation  amplitude  on  the  effective  damping  for  an  assumed  V-well 
representation  of  the  local  damping. 


SECTION  III 

RESULTS  AND  DISCUSSION 

The  polynomial  curve  fit  analysis  described  in  Section  II  has  been  applied  to 
experimental  nonlinear  pitch-damping  data  from  the  VKF  forced-oscillation  balance.  These 
data  and  the  acquisition  of  these  data  are  described  ip.  Refs  3  and  4.  Tabic  I  lists  the 
configurations  and  associated  test  conditions  that  were  considered  in  the  present  report. 
The  configurations  were  60-  and  70-dcg  half  angle,  blunted,  cones  with  various  afterbody 
arrangements  (Fig.  9). 

Preliminary  investigation  revealed  that  the  sixth-order  polynomial  representation  gave 
better  results  for  the  fitted  effective  damping  and  computed  local  damping  than  cither 
the  fourth-  or  eighth-order  polynomials.  All  the  results  presented  herein  are  from  a 
sixth-order  polynomial  curve  fit  unless  stated  otherwise.  The  polynomial  analysis  gave 
reasonable  results  for  the  majority  of  the  test  conditions  considered  for  these 
configurations. 

A  certain  amount  of  discretion  is  necessary  when  applying  this  technique  and 
evaluating  results.  A  good  polynomial  curve  fit  to  the  experimental  effective  damping  data 
does  not  guarantee  that  the  computed  local  damping  will  be  realistic.  That  is,  for  all 
practical  purposes,  the  method  does  not  yield  a  unique  solution  for  the  local  damping. 
A  slightly  different  fit  to  the  experimental  effective  damping  may  change  the  computed 
local  damping  considerably.  The  requirement  was  imposed  that  <he  computed  local  damping 
must  exhibit  realistic  behavior  before  the  results  can  be  considered  satisfactory. 
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TABLE  I 

VIKING  PITCH  DAMPING  TEST  CONDITIONS* 


Configuration 

M_ 

Red  x  10  6 

Figure  Number 

610N 

0.70 

0.570 

10a 

1.00 

0.440 

I  Ob 

1.30 

0  390 

10c 

1.50 

0  580 

10(1 

721 

0  70 

0.873 

11a 

0.80 

0.790 

lib 

0.90 

0.732 

lie 

1.00 

0.689 

lid 

1.10 

0.655 

lie 

1.20 

0.624 

Ilf 

1.40 

0.594 

i  is 

1.55 

0  581 

llh 

1.60 

0  579 

111 

1.90 

0.578 

Hj 

2.30 

0  624 

Ilk 

2  65 

1  299 

ne 

3.00 

0.730 

iim 

720 

0.70 

0  876 

!2a 

0.80 

0.791 

12b 

0  90 

0.729 

lie 

1.00 

0  686 

1 2d 

1. 10 

0  648 

12c 

1  20 

0  625 

12f 

1.40 

0  596 

l.’g 

1.53 

0  582 

12.' 

1.90 

0  587 

I2i 

2.30 

0  616 

i  2) 

3.00 

0  732 

1 2k 

*1  rom  Reft  3  and  4 
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60  deg  (Typ) 


5. 442  (Typ) 
l  -2.679 
x,^  *  1. 887 


1.450R  (Typ)  - 
Pivot  Axis  • 


10.000 
Diam  (Typ) 


Configuration  610N  =  601 


Xcg  *3.317 
l  -  2. 650  (Typ) 


Pivot  Axis 


70  deg  (Typ) 


3.875R  (Typ) 

All  Dimensions 
in  Inches 


15.500 
Diam  (Typ) 


Balance  Shield 


Configuration  720 


xcg  *3.317 


70  deg  (Typ) 


3.875R  (Typ) 


Pivot  Axis 


V  2.870  ; 

Hu  1 

/  11.388 

l - 1 

l  -2. 650  (Typ) 


15.500 
Diam  (Typ) 


Configuration  721 
Fig.  9  Model  Geometry 
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The  results  are  presented  in  Fig.  10  (Configuration  610N).  Fig.  11  (Configuration 
721),  and  Fig  12  (Configuration  720)  as  a  function  of  angle  of  attack.  Shown  are  the 
spread  in  the  experimental  data  (vertical  lines),  the  sixth-order  fit  to  the  experimental 
data  of  equivalent  damping  (uashed  line),  and  the  computed  local  value.  For  all  cases 
the  damping  is  stabilizing  at  higher  angles  and  varies  from  mildly  stabilizing  to  violently 
destabilizing  at  zero  angle  of  attack. 

Figure  13  presents  a  summary  of  the  damping  coefficient  at  zero  angle  of  attack 
versus  Mach  number.  The  basic  character  of  the  curve  is  the  same  for  each  configuration, 
i.c.,  near-neutral  st.  Mlity  at  a  Mach  number  of  0.7,  becoming  significantly  unstable  in 
the  Mach  number  range  from  1  0  to  2.2,  and  stable  at  higher  Mach  numbers 

Satisfactory  results  for  polynomial  representation  could  not  be  obtained  for 
Configurations  720  and  721  at  M»=  1.40  to  1.90  The  balance  output  for  these  conditions 
fFigs.  llg,  h,  i,  j,  and  12g,  h,  i)  reveals  that  the  experimental  effective  damping  has  an 
extremely  strong  (almost  discontinuous)  variation  between  o  of  zero  and  1  deg.  Also  shown 
in  these  figures  are  the  results  of  the  effective  damping  computed  from  an  assumed  V-wcll 
variation  of  the  local  damping.  Other  variations  of  the  V-well  and  square  well 
representations,  as  well  as  polynomial  curves,  were  also  investigated,  but  the  computed 
effective  damping  did  not  change  substantially  from  that  shown  in  these  figures  This 
computed  effective  damping  and  the  experimental  effective  damping  differ  considerably 
at  a  -  1  deg,  but  compare  favorably  at  a  =  0  and  2  deg 

For  M_  =  1.20  to  2.0,  it  is  possible  that  the  highly  nonlinear  behavior  observed 
in  the  experimental  effective  damping  could  be  attributed  to  a  destabilizing  hysteresis 
effect  in  the  static  pitching  moment  and  that  these  effects  are  greatest  at  a,.  =  0  (liq 
(29)).  The  mathematical  nature  of  the  curve  suggests  this  possibility,  but  there  is  no  relevant 
aerodynamic  theory  upon  which  to  base  this  speculation  Frce-oscillation  data  (for  various 
amplitudes  and  frequencies)  will  shed  some  light  on  this. 

SECTION  IV 

SUMMARY  AND  CONCLUSIONS 

A  procedure  to  extract  the  local  damping  from  the  experimental  effective  damping 
output  of  a  forced-oscillation  balance  has  been  derived  because  the  local  instantaneous 
pitch  damping  is  required  for  motion  simulation  of  certain  preliminary  Mars  unmanned 
reentry  (Viking)  configurations.  The  VKF  forced-oscili3tion  pitch-damping  balance  system 
has  been  briefly  reviewed  and  a  relationship  has  been  cstab!i<hcd  between  the  effective 
damping  output  of  the  balance  and  the  local  aerodynamic  damping.  Tins  relationship  is 
an  integral  equation  which  was  solved  to  yield  a  practical  connection  between  the  local 
and  effective  damping.  This  connection  involves  a  high-order  polynomial  (up  to  eighth 
order  in  a  and  )  representation  of  the  damping  coefficients.  Practical  use  of  this  technique 
requires  a  polynomial  curve  fit  to  the  experimental  effective  damping  which  results  in 
coefficients  for  the  polynomial  representation  of  the  local  damping.  This  procedure  has 
been  applied  to  experimental  data  obtained  by  the  VKF  forced-oscillation  balance  for 
preliminary  Viking  configurations. 
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The  following  characteristics  of  the  procedure  and  conclusions  concerning  its  practical 
usage  are  noteworthy: 

1.  Initial  attempts  to  obtain  a  practical  solution  for  the  integral  equation  (Eq. 

(31))  utilized  a  Fourier  series  representation  of  the  local  damping  and  also 
a  discrete  variable  approach.  Both  of  these  methods  were  unsatisfactory 
and  were  abandoned  in  favor  of  a  polynomial  representation  which  proved 
successful. 

2.  Discretion  is  necessary  when  applying  the  polynomial  curve-fitting  technique 
and  interpreting  the  results.  A  good  polynomial  curve-fit  to  the  experimental 
effective  data  docs  not  automatically  yield  reasonable  results  for  the 
computed  local  damping.  The  computed  local  damping  also  has  to  exhibit 
realistic  behavior  before  the  overall  results  can  be  considered  satisfactory. 

3.  The  sixth-order  polynomial  representation  gave  better  overall  results  for  the 
fitted  effective  damping  and  the  computed  local  damping  than  either  the 
fourth-  or  eighth-order  polynomials. 

4.  The  present  polynomial  technique  may  be  applied  to  frce-oscillation 
damping  data  acquired  by  oscillating  about  o  =  0.  if  do  is  set  equal  to 
zero  in  Eq  (39). 

5.  Tile  polynomial  method  may  also  be  applied  to  free-oscillation  data  obtained 
by  small  oscillation  about  a  nonzero  Oq.  In  this  case,  the  free-oscillation 
damping  must  be  determined  for  the  same  value  of  6j  from  the  oscillation 
envelopes  obtained  by  testing  at  different  values  of  a„. 

6.  If  the  static  moment  is  a  multivalued  function  of  a  for  an  oscillation  cycle 
(hysteresis),  this  will  create  additional  damping  (or  undamping)  which  will 
add  to  any  conventional  aerodynamic  damping  already  present  The  balance 
output  will  contain  the  sum  of  both  the  hysteresis  damping  and  conventional 
damping.  Forced-  or  free-oscilljtion  test  data  for  different  amplitudes  and 
frequencies  will  aid  in  interpreting  this  phenomenon. 

7.  The  apparent  damping  effects  of  a  multivalued  static  aerodynamic  pitching 
moment  are  not  included  in  the  integral  equation  (F.q.  (31))  relating  local 
and  effective  aerodynamic  damping.  This  also  applies  to  the  polynomial 
solution  of  the  integral  equation  (Eqs.  (34)  and  (39)). 

8.  The  polynomial  method  gave  reasonable  results  for  the  Mars  unmanned 
reentry  configurations  at  tl  •  isonic,  transonic,  and  supersonic  Mach 
numbers  above  M_  =  2.00.  .  *cvcr,  reasonable  results  could  not  be 
obtained  for  Configurations  720  and  721  at  Mach  numbers  from  1.40  to 
1.90. 
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a.  M_  =  0.70 


b.  M„  =  0.80 

Fig.  11  Local  and  Effective  Pitch-Damping  Coefficients  versus  Angie  of  Attack 
for  Configuration  721 
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C.  M_  =  0.90 


d.  M_  =  1.00 
Fig.  11  Continued 


Fig.  11  Continued 
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b.  M_  =  0.80 

Fig.  12  Local  and  Effective  Pitch-Damping  Coefficients  versus  Angle  of  Attack 
for  Configuration  720 


Fig.  12  Continued 
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